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0> . Abstract 
X 

I Tagged linearly polarized photons have been used at the Mainz Microtron 

I MAMI for simultaneous measurements of the p{^,p)'K^ and ^(7, 7r"'")n reaction 



channels to study the A(1232) transition. The energy dependence of 

the magnetic dipole Mj^^^ and electric quadrupole E^^^ amplitudes have been 
extracted from these data in the photon energy range from 270 to 420MeV. 
The E2/M1 ratio for the A(1232) transition has been determined to 

be —(2.5 ± O.lstat i 0.2sys)% at the resonance position (^33 = 90°). 
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I. INTRODUCTION 



Low energy electromagnetic properties of baryons, such as mass, charge radius, magnetic 
and quadrupole moments are important observables for any model of the nucleon structure. 
In various constituent-quark models a tensor force in the inter-quark hyperfine interaction, 
introduced first by de Rujula, Georgi and Glashow [^, leads to a rf-state admixture in the 
baryon ground state wave function. As a result the tensor force induces a small violation 
of the Becchi-Morpurgo selection rule 0, that the A(1232) excitation is a pure 

Ml (magnetic dipole) transition, by introducing a non-vanishing E2 (electric quadrupole) 
amplitude. For chiral quark models or in the Skyrmion picture of the nucleon, the main 
contribution to the E2 strength stems from tensor correlations between the pion cloud 
and the quark bag or meson exchange currents between the quarks. To observe a static 
deformation (ci-state admixture) a target with a spin of at least 3/2 (e.g. A matter) would 
be required. The only realistic alternative is to measure the transition E2 moment in the 
'jN A transition at resonance, or equivalently the E-^,^ partial wave amplitude in the 
A Nit decay. The amplitudes in the A^vr final state are usually denoted by E^^ and M^^, 
where E and M are the electric and magnetic multipoles, / is the isospin and / is the orbital 
angular momentum of the Nn system and the ± sign refers to its total angular momentum 
J = l± 1/2. 

The experimental quantity of interest to compare with the different nucleon models is the 
ratio Rem = E2/M1 = eI^_^ /mI^_^ of the electric quadrupole E2 to the magnetic dipole Ml 
amplitude in the region of the A(1232) resonance. In quark models with SU{Q) symmetry, for 
example the MIT bag model. Rem = ^ '^^ predicted. Depending on the size of the hyperfine 
interaction and the bag radius, broken SU{Q) symmetry leads to —2% < Rem < 
Larger negative values in the range —6% < Rem < —2.5% have been predicted by Skyrme 
models while results from chiral bag models give values in the range —2% to —3%. 
The first Lattice QCD result is Rem = (+3 ± 9)% and a quark model with exchange 
currents yields values of about —3.5% |jlO| . 
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The determination of the quadrupole strength E2 in the region of the A(1232) resonance 
has been the aim of a considerable number of experiments and theoretical activities in 
the last few years. Very recently, new experimental results have been published for the 
differential cross section and photon asymmetry of pion photoproduction off the proton 
from the Mainz Microtron MAMI and the laser backscattering facility LEGS at Brookhaven 
National Laboratory, with the results Rem = —(2.5 ± 0.2stat ± 0-'2sys)% from the Mainz 
group and Rem = —(3.0 ± 0.3stat+sys ± 0.2mod)% from the LEGS group ||12|. These new 



Rem results have started intense discussions about the correct way to extract the E2/M1 
ratio from the new experimental data. In particular the large variation in the Rem values 
obtained in theoretical analyses of these data at RPI [|l|] {Rem = -(3.2 ±0.25)%), VPI H 



{Rem = -(1-5 ± 0.5)%) and Mainz [jT5[ {Rem = —(2.5 ± 0.1)%) was quite unsatisfactory. 

In this paper we present the MAMI p(7, p)iT^ and p(7, TT~^)n differential cross sections and 
photon asymmetries and discuss in detail different analyses to extract the s- and p- partial 
wave amplitudes and the E2/M1 ratio. In Section II we briefly describe the experimental 
setup and show a selection of the measured cross sections and photon asymmetries for both 
pion production channels from the proton in Section III. The essential ingredients of a 
multipole analysis are outlined in Section IV. Three different analyses of our data to extract 
the s- and p-partial wave amplitudes and the E2/M1 ratio are discussed in Sections V and 
VI. We conclude with a summary and outlook in Section VII. 

II. EXPERIMENTAL SETUP AND DATA ANALYSIS 

Since the experimental setup used for this measurement was described in detail in Ref. 
p!6| , we will restrict the present discussion to the main features of the experiment. Linearly 



polarized photons were produced by coherent bremsstrahlung in a 100 fim thick diamond 
crystal [p!7| , p!8| . The photon energy was determined by the Glasgow tagging spectrometer at 
the Mainz Microtron MAMI |T^, which in a 352 channel focal plane detector analyses the 
momentum of the electron that has radiated the bremsstrahlung photon . This detector 
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system is able to energy tag photons in the range from 50 to SOOMeV with a resolution 
of about 2MeV The coUimation of the photon beam yielded a tagging efficiency of 

about 55% for incoherent bremsstrahlung. To continuously monitor the tagging efficiency 
and the photon polarization, a pair detector was used downstream of the hadron detector 
DAPHNE. This pair detector consists of a 0.5 mm thick Cu-converter followed by two 2 
mm thick plastic scintillators operated in coincidence. Its efficiency (e ~ 3% for photons) 
was regularly checked with a lead glass detector (e ~ 100% for photons) in calibration runs 
at low beam intensity. The photon polarization was determined from the photon spectrum 
measured by the tagging spectrometer in coincidence with the pair detector and with the aid 
of theoretical calculations [^]. The quality of these calculations was tested by an absolute 



measurement of the photon polarization using coherent 7r° photoproduction on ^He as a 
polarimeter reaction with an analyzing power A = 100%. Excellent agreement between 
calculations and experiment was found p3| and in this way, both the photon polarization 
and the photon flux could be determined with an absolute precision of better than ±2%. 

The liquid hydrogen cryogenic target was contained in a 43 mm diameter, 275 mm 
long Mylar cylinder with a wall thickness of 0.1 mm. The target density was stabilized and 
determined to an accuracy of ±0.5% by means of an an automatic pressure and temperature 
control system. 

The reaction products, the recoil proton from jp p7r° and the pion from jp nn'^, 
were detected using the large acceptance detector DAPHNE (21° < 9 < 159°, 0° < < 
360°) built by the CEA/DAPNIA-SPhN at Saclay and the INFN-sezione di Pavia |4 



Good definition of the charged particle tracks was obtained from the central vertex detector 
consisting of 3 coaxial cylindrical multiwire proportional chambers providing a polar angular 
resolution of A6' < 1° FWHM and an azimuthal resolution < 2° FWHM. This vertex 
detector is surrounded by a segmented — E — AE plastic scintillator telescope with 
successive thicknesses of 10 mm, 100 mm and 5 mm respectively. The outermost layer is a 
lead-aluminium scintillator sandwich designed to enhance the 7r° detection efficiency and to 
provide additional energy loss information for charged particles. 



In the first step of tlie analysis tliose events were selected, that had only one charged 
trajectory with its polar angle in the range 21^ < d < 159°. After this cut, the basic task of 
the data analysis was to identify the pions and protons. This separation was performed by 



using a range method as described in detail in Ref . , which simultaneously uses all of the 
measured energy losses in the scintillator layers of DAPHNE to identify charged hadrons 
(7r+, p) and to determine their energy. A restriction made on the vertex position, defined 
as the point on the reconstructed track, which hes closest to the detector axis, ensured a 
complete rejection of particles coming from the target windows and walls. As a result no 
empty target subtraction was needed. 

The main contributions to the systematic error in the determination of the unpolarized 
differential cross section are due to uncertainties in the photon flux (±2%), the target density 
(±0.5%) and (±2%) for the proton/pion separation with the range method. 



III. EXPERIMENTAL RESULTS 

In single pion photoproduction, the differential cross section for linearly polarized photons 
and unpolarized targets is given by the expression 

where S is the photon asymmetry and 6 and (j) are the polar and azimuthal angles of the pion 
with respect to the beam direction. Since DAPHNE has full 27r azimuthal coverage it allows 
a direct measurement of the (j) dependence of the differential cross section, and therefore, 
the determination of S and the unpolarized cross section dag/dQ at the same time. 

A selection of the measured cross sections and photon asymmetries in the cm-frame for 
the 7p piT^ and 7p —>■ rni^ reaction channels pT] , |26| are shown in Figs. This is the first 
data set for which the two observables dao/dQ and S have been measured simultaneously 
at all angles and photon energies for both pion production channels from the proton. The 
measurement covers the whole A (1232) resonance region in 16 energy bins between = 
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270MeV and 420MeV, where is the tagged photon energy in the lab-frame. The angular 
distributions for dao/dQ and S were measured from 45° up to 135'^ for the nvr^ channel. In 
the pTT^ case the angular range varied from 95° to 125° at = 270MeV to 65° to 125° for 
the highest energy point at E^ = 420MeV. In Figs. |1|-^ the result of a fixed-t dispersion 
relation analysis P^, which will be discussed below, is shown as a solid line. Our differential 



cross sections and asymmetries are in good agreement with the dispersion relation analysis 
except for the differential cross sections around 330MeV where we find a slight discrepancy. 

Figs. |] and ^ compare the energy dependence of the present MAMI results and the 
recent LEGS results [|l^ for 7p — > p7r° at 6^o = 90° and jp —>■ mc^ at 6^^+ = 85°. For 
the LEGS results the statistical error and the angle and energy dependent error have been 
evaluated point by point and combined to the net uncertainty bars, which are shown in 
Figs. ^ and ^ The error bars on the MAMI points are only statistical. The systematic 
scale uncertainties are ~ 3% for the MAMI results and ~ 1% for the LEGS results. There 
is an energy dependent discrepancy of up to 15% between the absolute differential cross 
sections. For both 7p p7r° and 7p nir^, the LEGS cross sections start to rise above the 
MAMI data at E^ = 280MeV, reaching ~ 15% higher values at the highest energy LEGS 
data point at E^ = 323MeV. Since the difference is energy dependent it will affect the 
absolute values of the resonant partial waves Mi/^ and (see discussion below). 

In Fig. 1^ the LEGS and MAMI angular distributions for the differential cross section 
and the photon asymmetry are compared at E^ = 320MeV for p{'i,p)TT^ and p(7, TT'^)n. One 
sees that the differential cross sections differ not only in absolute magnitude, but also in 
angular shape. The two most forward angular points at Ot^o = 70° and 80° of the LEGS 
]97r° differential cross section drop faster than the MAMI results. The same behavior is seen 
for the n7c~^ cross section at 6*^+ = 20°, 150° and 170°. This difference in the shape of the 
differential cross sections for the two data sets will play an important role in the discussion 
of the non-Born contribution for partial waves with > 2. On the other hand, the photon 
asymmetry agrees well for the pir^ channel and shows only small differences for the nvr"*" 
data, where the MAMI results are slightly above the LEGS data. 
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IV. MULTIPOLE ANALYSIS OF THE MAMI DATA 
A. General Aspects 

In 'fN Nit both incident particles and the final nucleon have two spin states yielding 
in eight degrees of freedom. From parity conservation this is reduced to a total number of 
four complex amplitudes to describe the reaction. Allowing for one arbitrary phase factor, 
there are therefore seven independent physical quantities which need to be measured at any 
setting of and 6-,^. Different sets of amplitudes are used in the literature, e.g. the helicity 
amplitudes {Hi, H2, H3, H4) introduced by Jacob and Wick or the Pauh amplitudes 
(Fi, F2, F3, F4) by Chew, Goldberger, Low and Nambu |2D|. While the observables of single 
pion photoproduction are more elegantly expressed by the helicity amplitudes, the Pauli 
amplitudes are particularly suited for a decomposition into partial waves |]30| , |3l| . 

A complete database for pion photoproduction (the "complete" experiment) requires at 
least 8 independent observables to specify the multipole amplitudes to all orders in |]3^-|34|. 
Such complete information is not available at present and standard multipole analyses have 
to rely on the differential cross section da/dQ and the three single polarization observables 
S (photon asymmetry), P (recoil nucleon polarization) and T (target asymmetry). For pion 
photoproduction from threshold up to the A(1232) resonance region, these four observables 
provide sufficient conditions for a complete database, if the higher partial waves /jr > 2 
can be adequately represented by the Born contributions. This approach is expected to be 
appropriate up to the region of the higher resonances [E^ ^ 600MeV) and even there may 
affect only certain multipoles (e.g. -Di3(1520) ^ E2~, M2-). Such arguments were used 
by V. Grushin to analyse the Kharkov data |^ {da/dQ, S, P and T) for ■yp pn^ 
and 7p —>■ rni^. In their analysis both the real and imaginary parts of the s and p wave 
amplitudes could be determined for the first time independently of the pion-nucleon phase 
shifts . The mean difference between the 533^^^ phase from the Kharkov analysis and 
the known pion-nucleon scattering phase 5"^^ was only —(2.3 ± 0.5)° over the entire energy 
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range = 250 - SOOMeV. 

The pn^ and mi~^ multipole amplitudes A4i± {Mi± stands for Eq+, Mi-, . . .) are related 
to the isospin 1/2 and 3/2 components and Ai^^ by 

M.^ipn') = M]i' + Im%' (2) 
Mi^inn^) = V2{M]i' - lu'ii') . (3) 

The multipole amplitudes Aij^ are complex functions of the cm. energy W. Below the two 
pion production threshold {E^ ~ 310MeV), the Fermi- Watson theorem [^] allows one to 
express the phases of the complex multipole amplitudes by the corresponding pion-nucleon 
scattering phase shifts 

Ml^ = \MUexp{t5l^ + nn), (4) 

where / can be 1/2 and 3/2 and n is an integer. However, even at energies of about 400MeV, 
the irN inelasticities in the P33 partial wave are very small, which suggests that the Fermi- 
Watson theorem can be applied well above the two pion threshold. 

There exist two basic methods to extract the multipole amplitudes from the database, 
the "energy independent" and the "energy dependent" approach. In the energy independent 
approach, each energy is investigated independently by use of standard minimization 
techniques, the fit parameters being the real and imaginary parts of the multipole amplitudes 
Adjj.. Below the two pion production threshold the Fermi- Watson theorem is used, which 
reduces the number of the necessary observables for a complete experiment by a factor 
of two, because only the absolute values |A^h-| of the partial wave amplitudes need to 
be determined from the fit. In the energy dependent approach, the data at all energies 
are analyzed simultaneously. Either an energy dependent parametrization of the partial 
wave amplitudes must be assumed, or as in our case, the energy dependence is taken from 
dispersion relations. The principal advantage of this method is that continuity is built in 
from the beginning and systematic errors tend to cancel out. In some cases one or the other 
of these two approaches is to be preferred. If the data are closely spaced in energy and 
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cover both differential cross sections and polarization observables ("complete experiment") 
at each energy, the energy independent approach is more advantageous. If on the other 
hand the data are widely spaced and only few polarization observables are available, the 
energy dependent approach is the better one. This approach is also useful if the general 
resonance structure is already known and the main interest is to obtain the small partial 
wave amplitudes. However, the experimental data has to cover the complete energy range 
of the dominating resonances in order to allow a reliable separation. For instance in the 
case of the Mi+ amplitude in the — A(1232) transition, one has to cover the complete 
resonance region (250 — 450MeV) to get the Mi+ multipole as precisely as possible l3^j3^ . 



B. Higher partial waves 

A general problem common to both the energy dependent and the energy independent 
approach is to decide at which value of the angular momentum 1^^ the partial wave expansion 
should be truncated in the fit and how to approximate the higher partial waves. The 
maximum value rimax that the data can determine is found by fitting the angular distributions 
to a Legendre series or, equivalently, a power series expansion in cos6 



t'TTtaX 



S^-^ = fE ^Sco."^. (6) 
ail sm^B k 

r^_i^ = f V (7) 

dVLsinO k ^ ^ ^ 



n- 



max 



P^-L- = i V A^cos^^e. (8) 

The experimentally accessible polynomial coefficients A" are quadratic or bilinear products 
of the electric Ei± and magnetic Mi± multipole amplitudes. 

In this section we examine the sensitivity to higher partial wave contributions of the 
coefficients A" and extracted from the present data. The difference of the cross sections 
for the photon polarization perpendicular and parallel to the reaction plane, da±/dQ — 
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da\\/dQ, is particularly sensitive to higher partial waves. With partial waves up to = 2 
the difference 

1 f da± da\\\ ^da 1 q.^ n ^ in\ /^\ 

S— — — = f-(As + 5sC0S^ + CEC0s2^) (9) 



2sin^e \dVL dVt ) dVt sin^O k 
has three polynomial coefficients A-^, By; and Cs with 

— ^Y:{Swavey Pwave) + R-G [-E'o+'^tya»;e] + M«)ai)eP ; (lO) 

i?s^Re[(Mi+-Mi_)rf:_] , (11) 

wave 

(12) 

In the case that only the s wave (£"0+) and p waves (Mi+, Mi_ and £"1+) contribute, this 
difference would be proportional to and therefore independent of the pion angle 6. The 
multipole dependence of the coefficients A^, and is described in Appendix A. 

Fig. |] shows the difference of da±/dQ — da\\/dQ for both pion reaction channels on the 
proton at our lowest photon energy = 270MeV, at = 350MeV and at our highest 
energy E^ = 420MeV. In the pir^ channel one recognizes only above the resonance a small 
deviation from the constant behavior. This is due entirely to the Born contribution to the 
By coefficient, i.e. the interference between the real part of the dominant Mi+ amplitude 
and the real part of the d wave amplitudes, for example in terms like Re(Mi_|. — Mi_)Rei?2-- 
Such contributions become extremely small at resonance since Re(Mi+ — Mi_) ~ 0. 

The behavior of the difference, da±/dfl — da\\/dQ, for the 7p —>■ pn^ data from the LEGS 
collaboration is different. Their observed angular dependence seems to require a sizeable 
non-Born contribution from d and / waves. The origin of the different angular dependence 
of the LEGS results for da±/dQ — da\\/dfl arises from the different shape of their unpolarized 
cross section data. The LEGS photon asymmetries alone are well described by our multipole 
fit which takes Born contributions in higher partial waves of > 2 into account. 

In contrast the difference of da±/dQ — da\\/dVt for 7J9 rm^ shows a strong angular 
dependence at all energies. The main reason is that in charged pion photoproduction the 
pion pole Born graph leads to significant contributions of higher partial waves. 
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The sensitivity of the differential cross section to higher partial waves is most pronounced 
at the extreme forward and backward angles. Taking partial waves up to = 2, the 
differential cross section is 

^ = |(A + Bcos{e) + Ccos\e) + Dcos^ie) + Ecos^ie)) (13) 

with five coefficients of the form 

^ — -^(^Siiiave^ Pwave) ~l~ [.Eo+C^ujai'e] ~l~ |'J^uiat)e| j (-^^) 
B ~ B{s^ave,Pn^ave) + Re [(Mi+ - M,_)dl,J , (15) 
C — C ^Swavei Pwave) ~\~ [-Eg+ciuiarie] ~l~ |c^«)atie| ) ('^^) 

D~Re[(Mi+-Mi_)rf;,,e] , (17) 

E = \d^i,ave\'^ ■ (18) 

The effect of the d waves is largest for the coefficients B and D, because of an interference 
term between the large Mi+ amplitude and the d waves. However, the contributions of 
these terms can be neglected at the top of the resonance (^33 = 90*^). As an example this is 
illustrated for the combination 



Re 



(Mi+ - Mi^)E*_ = Re(Mi+ - Mi_)Re£'2_ + Im(Mi+ - Mi_)ImE2- , (19) 



where the first term vanishes, because Re(Mi+ — Mi_) passes through zero at the resonance 
position (^33 = 90°) and the second term can be neglected, because Imi?2- is small due to 
the small d-wave phase. Fig. ^ illustrates the sensitivity of the differential cross section to 
higher partial waves. The solid line shows the ratio {da'^'^ / da^^^^) of the differential cross 
section for only s- and p-wave contributions to the differential cross section with higher 
partial waves. The calculated cross section ratio is shown at the pion angles O^^o = 0^, 90° 
and 180° in the energy region between 200 and SOOMeV. For this calculation ||39|| the SAID 
solution SM95 was used, in which partial waves up to = 3 were allowed to vary to fit 
the data set but the normal Born contribution was assumed for l.^ > 4. At ^^0 = 90°, the 
contributions from the higher partial waves are far below 1%, because they arise only from 
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an interference term with the small s wave Eq+, for example in the form Re(i?o+'i^a„e) as in 
eq. 15. Below and above the A(1232) resonance, however, contributions from /^r > 2 are of 
the order of 10 — 20% of the differential cross section at 0° and 180°. These contributions 
arise almost completely from the interference between the dominant Mi+ amplitude and the 
Born parts of the higher partial waves. In our analysis there are no indications of significant 
non-Born contributions for higher partial waves (/tt > 2). 



C. Influence of systematic errors in the database 

An important issue for any multipole analysis is the question of how to handle the sys- 
tematic errors of the experimental results for the different observables used in the database. 
The ideal database should contain experimental results for both 7p — pn^ and 7p nTr"*", 
measured simultaneously with one setup for all photon energies and covering the full range 
of polar angles to minimize energy and angular dependent systematic errors. The analysis 
of such data would yield a consistent separation of the isospin 1/2 and 3/2 parts from the 
partial waves of the pn^ and ravr"^ reaction channels. However, the isospin separation will fail 
if there is a considerable inconsistency, such as a normalization error, between the 7p mr^ 
and 7p pn^ observables {da/dQ, S, T, . . . ). The existing pion photoproduction data, for 
example the SAID database are far from the requirements of the ideal case. Most of the 
experimental data are taken at one photon energy at a time, the angular distributions are 
not measured simultaneously for all angles and, more importantly, the consistency between 
different measurements is less good than their assigned errors would suggest. When combin- 
ing data from different experiments into one database, one has the difficult task to account 
properly for the angular and energy dependent systematic errors. A common procedure is 
to combine the experimental systematic errors in quadrature with the statistical errors or 
to multiply all data from one set with a systematic scale error {cgys) by a common factor 



/ while adding (/ — 1) /a^y^ to the x E3- '^^^ latter method allows for adjustments of 



an overall scale (angle independent systematic errors), but not for systematic effects in the 
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shape of the angular distributions which is important for the determination of the small E2 
amplitude. As shown in the following paragraph, it is the shape of the differential cross 
section and the photon asymmetry near 6 = 90^, however, which is sensitive to the E2 
amplitude. Combining all existing pir^ and rnr^ data in one database regardless of their 
consistency will result in mean values for the dominant partial waves (Mi+ and -Eo+) but 
meaningless values for the small partial waves and Mi_). Any reliable partial wave 

analysis for the small amplitudes will require a careful selection of the experimental results 
used in the database. In practice one has to choose observables that are sensitive to the small 
amplitude from one experiment, or one has to use analyzing methods which are insensitive 
to energy and angle dependent systematic uncertainties. 



Fig. |T0| shows da /dVL and S for different multipole contributions at the A(1232) resonance 
((^33 = 90°). This illustrates the sensitivity to the higher partial waves (/jr > 2) and to 
the small amplitude. The contribution of the higher partial waves to the 7p mT~^ 
channel is significant at all angles. For the differential cross section at 6^^+ = O" and the 
photon asymmetry around = 90° the full calculation and the truncation to s and p 
waves differ by as much as a factor of two. The 7p — pn^ case is different and one cannot 
distinguish between the full calculation and the result where only s and p waves are taken 
into account. The main reason is that the pion pole term does not contribute in neutral 



pion photopro duct ion. Fig. |T0| shows also the sensitivity of S and da / dVt to the small Eij_ 
amplitude for both pion channels. The main effects of the Eij^ multipole in the pvr^ channel is 
found near ^^o = 90° for the photon asymmetry and near 0°, 90° and 180° for the differential 
cross section. This sensitivity has its origin in the interference term Re(£'i+M*_,_), which 
appears in the coefficients 

A ^ ^|Mi+|2 - 3Re(Ei+M*+) , (20) 
C ~ -|mi+|2 + 9Re(Ei+A'/i%) , (21) 



2' 

3,,. ,2 



|Mn.|' + 3Re(Ei+M*+) , (22) 
and contributes to both da/dQ and S. Without Ei^ strength the shape of the differential 
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cross section would be determined by C/A = —3/5; and the photon asymmetry by S(90°) = 
Ay.IA = 3/5. For E1+/M1+ = -2.5% these values would change to C/A = -0.7 and 
S(90°) = 0.54. Setting the multipole to zero (dotted line) results in an increased 
differential cross section at the extreme angles O" and 180° and a lower value at 90*^. The 
photon asymmetry in contrast is enhanced at 90° for Ei^ = 0. The same behavior for S is 
found in the nvr"*" channel. In the differential cross section for rnr^ the effect of EiJ^ is most 
sizeable at the backward angles. 

From the above discussion it becomes obvious that a reliable extraction of the small £'1+ 
amplitude will require precise photon asymmetry data around 90° and precise data for the 
differential cross section at all angles. 



V. DISCUSSION OF THE DIFFERENT ANALYSES 

We have performed three different analyses of our data in order to extract the s- and p- 
wave amplitudes and the E2/M1 ratio. First, a fit to the differential cross section and photon 
asymmetry only in the '~fp pir^ channel; second, an energy independent simultaneous 
multipole analysis of the '-/p pn'^ and 7p — ^ nn^ data to get the isospin separation in 
the whole energy region (270 — 420MeV). In the third analysis an energy dependent fixed-t 
dispersion analysis has been performed, which includes additional observables in order to 
study the stability of the different multipole solutions. 



A. Analysis of the MAMI pvr" data 

Since the result from the fit to the 7p pn^ channel alone has already been pub- 
lished [|ri|], we summarize here only the key points of our analysis. The p7r° angular distri- 
butions for the unpolarized cross section dao/dQ, the parallel part da\\/dVL (pion detected 
in the plane defined by the photon polarization and the photon momentum vector), and 
perpendicular part da±/dQ can be expressed in the s- and p-wave approximation by the 
parametrization 
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= + B,cos{d) + C,cos\e)) , (23) 

where q and k denote the center of mass momenta of the pion and the photon, respectively, 
and j indicates the parallel (||), perpendicular (±) and unpolarized (0) components. The 
coefficients Aj, Bj and Cj are quadratic or bilinear functions of the s- and p-wave ampli- 
tudes. In particular, da\\/dQ is sensitive to the Ei+ amplitude, because of interference with 
Mi+ in the terms 

A|| = I Eo+ P + I 3Ei+ - Mi+ + Afi_ |2 , (24) 
fill = 2Re[Eo+(3^i+ + Mi+ - Mi_)*] , (25) 
Cii = 12Re[Ei+(Mi+ - Mi_)*] . (26) 

Furthermore, the ratio 

o ^ 1^ ^ Re(Ei+(Mi+-Mi_)*) 

12 A|| I Eo+ |2 + I 3^1+ + Mi+ - Mi_ |2 ^ 

is very close to the ratio Rem = /^i+ as we will show below. At the A(1232) resonance 
{^33 = 90°) the ratio R can be identified with the Rem value. 



ImEi(^ 

R — Rem — TjK 

ImM^y 



(28) 

W=Ma 



This is the key point of our analysis. The method has the advantage of being independent of 
absolute normalization and insensitive to many systematic errors, because Rem is extracted 
from the ratio of the coefficients C\\ and A\\ fitted to the angular distribution of da\\/dQ. 



Analyzing the data by this method the following result is obtained |TT 



Rem = -{2.5 ±0.2,tat)%. (29) 

Without any correction for neglecting the isospin 1/2 contributions to Ei^ and Mi+, the 
maximum absolute systematic error for Rem would be ±0.5%. To further reduce this 
systematic error one can take the Born contribution and estimate the size of the isospin 1/2 
contribution for £'1+ and Mi+. With a conservative estimate the absolute systematic error 
reduces to less than ±0.2% (see Appendix B for details). 
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B. Energy independent multipole analysis of the MAMI pTr^ and nvr"'" data 

To obtain the isospin decomposition for the s and p waves in the whole energy range 
(270 — 420MeV), we have performed an energy independent multipole analysis of our data. 
In this analysis we fit 8 parameters, the absolute values of the s and p waves for isospin 1/2 
and 3/2 (Eo(^ E^i^ mI'^ , mI'^ , mI'^ and mI'^), to the photon asymmetry 

and the differential cross section for the 7p pn^ and 'jp nn^ reaction channels. Higher 
partial waves (Zyr > 2) are taken into account by the Born terms, including p and u exchange 
in the t-channel. We obtain the coefficients A, B and C from the angular distribution of 
the differential cross section and A-^ from the photon asymmetry. The two pion reaction 
channels are decribed by 8 coefficients, which are independent combinations of the isospin 
1/2 and 3/2 multipole amplitudes of the s and p waves. In addition the Fermi- Watson 
theorem is used to determine the real and imaginary parts of the s- and p-wme multipole 
amplitudes. 

In Fig. [Tl| we show our fitted multipole amplitudes as ReA^fj_/cos5fj_, which is the abso- 
lute value of the amplitude up to a sign that can change with photon energy as in the case of 
Ei^_^. The values of the s- and p-wave amplitudes are fitted independently at each energy to 
the MAMI data 7p pn^ and 7p —>■ nir^. The results for the isospin components for the s- 



and p-wave amplitudes are presented for 16 photon energy bins |]2B[. Although systematic 
errors are greatly reduced because the pvr^ and mr^ data are obtained with the same set up, 
this method is sensitive to small remaining systematic differences in the results of da/dQ 
and S for the two pion channels, because it relates 8 parameters to the 8 "observables" 
(A,B,C and A^ coefficients for the two pion channels). Including additional polarization ob- 
servables in the database and performing an energy dependent multipole analysis overcomes 
these shortcomings. 
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C. Energy dependent multipole analysis of the MAMI data 



In a third analysis the MAMI pion photoproduction data are analyzed using fixed-t 
dispersion relations based on Lorentz invariance, isospin symmetry, unitarity and crossing 
symmetry [|15],^. This analysis includes the recent MAMI data for the differential cross 
section da/dQ and photon asymmetry S for pn^ and mr^ from the proton p6| , ^ , |47[| , both 



older and more recent data from Bonn for the target asymmetry T ||48|-^ , and differential 



cross section data on vr production off the neutron from Frascati and recently from 



TRIUMF 52 



With this method we performed both an energy dependent and an energy independent 
multipole fit as shown in Fig. |1^ by open circles (energy independent) and the solid line 
(energy dependent). The agreement between the E^^^ and Mfi^ multipole results obtained 
from the dispersion analysis and the multipole analysis described in the previous section 
is very good. Both energy independent results show the same energy behavior for the two 
multipoles. For the real and imaginary parts of the amplitude the two results are 

in excellent agreement, the solid dots lying on top of the open circles up to 425MeV (see 
Fig. |l^). The agreement is also very good in the case of the small E'^!^^ amplitude. In 
particular, we find that the inclusion of other observables in addition to the MAMI pion 

3 /2 

photoproduction data decreases the fluctuations slightly for the E^+ multipole. 



The results of the energy dependent fit (solid line in Fig. |T^) agree well with the Mi+ 
multipole derived from the energy independent fit, but small systematic differences may be 
seen for the electric multipole at energies above the A(1232) resonance {E^ > 360MeV). 
The ratio Rem at the resonance agrees quite well in both dispersion analyses. 

Rem = -{2.5 ±O.Uat)% (30) 

for the energy dependent fit and 

Rem = -{2.33 ± 0.17 stat)% (31) 

for the energy independent fit. 
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In addition we have checked the effect of changes in the database. Replacing our differ- 



ential cross section by the Bonn data obtained in the seventies [0,0] reproduces the results 
for the leading multipoles but changes the smaller multipoles to some degree (see Tab. |). 
The Bonn results for da/dQ at very forward and backward angles force the fit to a smaller 
Ei^ value near the resonance. However, it should be kept in mind that the overall compi- 
lation of the Bonn 'yp pix^ data by Genzel et al. |Q results from different experimental 



setups, which have been re-appraised and re-evaluated before combining them into one data 
set. The main part of the angular distribution (50° < 6^ < 160°) for dcr/dQ was measured 



by detecting the recoil proton in a magnetic spectrometer |55|, i.e. by fixed angle and single 
energy measurements. The differential cross section for 10° < < 70° was measured by 
detecting the two decay photons of the 7r° with lead glass blocks in a different experiment 
||56|| . Each of these experiments has a systematic error of the order of ±5%, and it takes 
much less than this systematic error to change the shape of da/dQ {C/A ratio) and with 
this the Rem value. 

In the case of the LEGS data the larger cross section leads to larger values for the 
resonant multipoles Mi!^ and E'i^^. The main reason for this is that the cross section 
difference between the LEGS and MAMI data is energy dependent (resonance behavior) as 
shown in Fig. ^ and Fig. |^. Inclusion of the polarization data T, S and P from Kharkov 
pGl does not affect our fit because of the large statistical and systematic errors especially 
for the 7p — > p7r° channel. 

VI. DISCUSSION OF THE RESULTS 

According to the Fermi- Watson theorem the and M^^ partial waves have the same 
phase and the ratio Ei^^/M^!^ is a real quantity. As shown in Fig. 0, this ratio is 



strongly dependent on the photon energy and varies from —8% at E^ = 270MeV to +2% 
at E^ = 420MeV. The ratio Rem is defined at the A(1232) resonance position, where 
533 = 90°, by 
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p3/2 

^EM — TTT 



(32) 



We note that within the K-matrix fomahsm this ratio is free of background contributions 
since is the K-matrix pole. The extraction of the genuine A(1232)-resonance parts of 
the magnetic dipole and electric quadrupole multipoles is model dependent, see for example 
the work on dynamical models [p8H60(| and the recently proposed speed plot analysis ||T5| , |6l| . 
In order to get a reliable value for Rem it is important to use a database with small angle 
dependent systematic errors, because the small amplitude depends on the shape of 
da/dQ and the absolute magnitude of the photon asymmetry at 90°. To obtain a consistent 
separation for the isospin 1/2 and 3/2 parts of the partial waves, experimental data are 
required with small inconsistencies in the absolute normalization between the 7p pn^ and 
7P — > nn^ observables. To reduce the influence of such systematics in the multipole analysis, 
the observables S and da/dQ for both reaction channels pn^ and mr^ were simultaneously 
measured at all angles and energies with the DAPHNE detector at MAMI. The three different 
analyses to extract Rem from the MAMI data agree very well with each other, and at the 
A(1232) resonance we get the final result 

„ Imi?i ( , . ^ , , 

Rem = = - 2.5 ± O.Uat ± 0.2,^,, % . 33 

In addition, the discussions about the extraction of the correct Rem ratio from our data, 
which arose in the literature after our first publication |jll| can now be summarized as follows. 



In the Comment of the VPI group [14], our pre data were included in the SAID database and 



a value Rem = —(1.5 ±0.5)% was obtained. As pointed out in our Reply ||41[] the difference 



between the VPI result and our value Rem = —(2.5 ± 0.2)% is due to the database used 



in the SAID analysis. As has been recently confirmed by the VPI group |^2|, the exclusion 
of all pre-1980 pn^ differential cross section data in the SAID database changes Rem to 
—2.5%. It is obvious that the key problem arises from the inconsistencies in the different 
data sets used in the analysis. The forward and backward angles for the differential cross 
sections used in the SAID database force the VPI solution to smaller Rem values. However, 
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this solution overpredicts the Mainz photon asymmetry data at resonance, S(90°) = 0.58, 
compared to the Mainz result, S(90°) = 0.54. The reason is that at resonance (^33 = 90°) 
the fit has only one parameter, the Rem value, to describe simultaneously the shape of the 
differential cross section 

C -3/2 + 9Rem 



A 5/2 - 3Rem 

and the photon asymmetry at 90° 

Ay, 3/2 + 3Rem 



(34) 



(35) 



A 5/2 - 3Rem ■ 

As we have stated before, the compilation of all the existing experimental data into one 
database can result in mean values for the dominant multipoles and meaningless values 
for the small multipoles. For example, if one combines the new MAMI data with the 
Bonn differential cross section results, then the Rem value is affected by small systematic 
differences in the two data sets. The Bonn differential cross sections range from 10° to 160° 
while the MAMI cross sections cover angles from 75° to 125° and are slightly below the Bonn 
data at = 340MeV. In the combined data set the shape of the differential cross section 
is changed, i.e. the C/A ratio gets smaller. This is the main reason that for the combined 
Bonn and MAMI data set the Rem value is even below the value which one gets from the 
Bonn data alone. This influence of the Bonn data on the extracted value of the ratio Rem 
has been confirmed by the BNL and RPI/VPI groups . 



In a second Comment [|1^ our ]97r° data were analyzed by the RPI group, who obtained 
the result Rem = —(3.2 ± 0.25)%. However, the inclusion of our mc^ data in the database 
lowered this value to Rem = —(2.64 ± 0.25)% [^, in agreement with our analysis. 



VII. SUMMARY AND OUTLOOK 



We made the first simultaneous and accurate measurements of the differential cross sec- 
tions and photon asymmetries for the reactions p(7,p)vr° and p{'j,7T~^)n in the A(1232) 
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region. Using this data we have performed a multipole analysis to obtain the isospin de- 
composition of the s- and p-wave multipoles and extract the E2/M1 ratio over the energy 
range 270 — 420MeV. Our final results at the A(1232) resonance position (^33 = 90°) 
are Rem = -(2.5 ± 0.1,^* ± 0.2,^,,)%, and Ay2 = -(131 ± l){10-yVG^) and A3/2 = 
-(251 ± l){10-^/VGeV) for the helicity amplitudes. 

In the meantime more ^(7, TT^)p data for the differential cross section and the photon 
asymmetry have been taken. This new experiment covers the full range of polar angles by 
observing the two decay photons in the TAPS detector at MAMI. The analysis is in progress 
and will produce new differential cross sections for the extreme forward and backward angles. 
In this way we hope also to clarify the absolute normalization problem in the differential 
cross section between LEGS and MAMI. 
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Appendix A 

The coefficients A, B, C and A-^ in the s- and p-wave approximation for the differential 
cross section and the photon asymmetry: 

A{s^ave:Pu,ave) = l^o+l' + ^l^i+l' + ^iMi+f - \Mi_\^ - 3Re(Ei+(Mi+ - Mi_)*) + Re{M,_M^^) 
B{s^ave,Pwave) = 2Re{Eo+{3Ei+ + Mi+ - Mi_)*) 

C{s^ave,Pv.ave) = ^l^l+l' " ^l^i+j^ + 9Re(Ei+(Mi+ - Mi_)*) - 3Re(Mi_M*+) 

As{s^ave,Pwave) = -^l^^i+l' + + 3Re(£;i+(Mi+ - Mi_)*) + 3Re(Mi_Mi%) 



The differential cross section up to /^r = 2 partial waves 

^ = ^(A + Bcos(0) + Ccos^e) + Dcos\e) + EcosHO)) 
d\l k 

has five coefficients with 

A — A(^Sy]a,veiPwave) 

5 45 g g 

+ -\E2-\^ + Yl^2+r + -|M2_|' + -IM2+P + Re(Eo+(-6S2+ - 3M2+ - + 3M2_)*) 

+Re(£;2_(y£;2+ + 3M2- + 3M2+)*) + Re(M2_(-9M2+ - ^-E2+y) + Re(^M2+£;2%) 

B — BiyS-ujaveiPwave) 

+2Re(Mi+(£;2_ - %E2+ + 9M2+ + 6M2_)*) + Re(Mi_(18£;2+ - 2£;2_ + 9M2+ + 6M2_)*) 
+6Re(^i+(-2E2- + 3M2_ + 3M2+ + 3^2+)*) 

C — C (^Syjavei Pwave) 

-^1^2-1' + y 1^2+1' + ^|M2_P + 27IM2+P + Re(Eo+(6E2+ + 3M2+ + ^2- - 3M2_)*) 
+Re(£;2-(-63£;2+ - 9M2_ + 9M2+)*) + Re(M2_(81£;2+ + 8IM2+)*) - Re(81M2+£;2%) 

D = 2Re(Mi+(30E2+ - I2M2+ - 18M2_)*) + Rc(Mi_(-30E2+ - I5M2+)*) 

+Re(£;i+(-36M2- + 36M2+ + 18£;2+ + 18£;2-)*) 

45 45 45 

E = — |^2+r - y|M2+|2 - 90Re(M2+M*_) + yRe(E2+(5M2+ - bM^- + SE^-)*) 
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The photon asymmetry E 
has three coefficients with 

— ■^'E^Syjavei Pwave) 

-^\E2-f - 9|^2+r - l\M2-\^ + Re{Eo+{E2- + E2+ + M2_)*) 
+Re(£;2_(3M2_ + '^E2+)*) - ^Re(£;2+M*_) 
5s = +3Re(Mi+(5E2+ + 4M2+ + 6M2_)*) - 3Re(Mi_(5E2+ - 5M2+)*) 

-3Re(£;i+(9£;2+-9£;2-)*) 

45 45 
= — IM2+I' - 45|E2+r + — Re(E2+(3E2_ - 2M2- - M2+)*) + 90Re(M2+3M*_) 
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Appendix B 

At the A(1232) resonance position, where the phase ^33 passes through 90° {E. 
340MeV) we find ReMi+(3/2) = 0, Re(Mi+ - Mi-) ~ and neghgible contributions from 
higher partial waves {l-,^ > 2) (see sect. IV b). The ratio R of the coefficients A\\ and C\\ for 
p(7,p)7r° can then be expressed by 

12A\\ |Eo+P + |3^i+ - Mi+ + Mi_|2 ■ 

Neglecting l-Bo+P, 9|£'i+|^ and all terms with Re(Mi+ — Mi_) results in 

^ _ ImSi+Im(Mi+ - Mi_) 

" -6Im^i+Im(Mi+ - Mi_) + lm^{Mi+ - Mi_) ' 

Dividing by Im^(Mi+ — Mi_) we obtain 

1 — 6R^o ' 

with 

IiiliTi- 



ImMi+ - ImMi_ 

If we further neglect ImMi_ and the isospin 1/2 components ImM^^^ and lmEl!^^, then the 

1 /2 

largest correction comes from ImEi!^ , which is of the order of 10 — 20% of Im£^i+. The final 
result depends on the ratio of Im£^i{.^/Im£'i+(p7r°). If we further assume it! = —2.5% we 

obtain: 



(a) for ImEl!^ = 



C|| I-ReM r, o nco/ 



(b) for Im^i(Vlm^i+ = 10% 



^=T^" /"Rt^" ^ i?.M = -2.65%, 
12^11 1 — o.oRem 



c) for ImEi{.7lm^i+ = 20% 
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with 

Without any correction for the isospin 1/2 contributions to Ei^ and Mi+ the maximum 
systematic error of Rem would be ±0.5% absolute. To further reduce this systematic error we 
can estimate the size of the isospin 1/2 contribution by the Born contribution. A conservative 
estimate of the remaining systematic error is less than ±0.2% absolute. 
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TABLES 



Data 



ImM' 



3/2 



Im£;: 



,3/2 



Rem{%) 



MAMI 
LEGS 
Bonn 



37.66 
39.35 
38.31 



-0.924 
-L276 
-0.643 



-2.54 
-3.24 
-L68 



TABLE L Imaginary parts of the multipoles E^^_^ and M'^^_^ in units of 10 '^/m^+, and the 



ratio Rem = ^"^E^^-^ /lmM'l'^\w=Mi^ in the energy dependent dispersion analysis. The different 



f3/2 



rows correspond to differential cross section data from MAMI |26,4^, LEGS |12| and Bonn |53,54|. 
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FIGURES 
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FIG. 1. The differential cross section da/d^l for p(7,p)7r° at 6 different photon energies. The 
soUd hne shows the result of the energy dependent multipole analysis. 
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FIG. 2. The photon asymmetry S for p{'y,p)TT^ at 6 different photon energies. The sohd Une 
shows the result of the energy dependent multipole analysis. 
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FIG. 3. The differential cross section da/dVt for p(7,7r^)n at 6 different photon energies. The 
solid line shows the result of the energy dependent multipole analysis. 
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FIG. 4. The photon asymmetry S for p{^,iT^)n at 6 different photon energies. The sohd hne 
shows the result of the energy dependent multipole analysis. 
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FIG. 5. The energy dependence of the differential cross section at 9j^o = 90° for p{'y,p)Tr^ in 
comparison with our energy dependent multipole analysis (full line) and the VPI solution SP97k 
(dotted line). For the LEGS points the statistical error and part of the systematic error are 
combined to the net uncertainty bars |^]. The error bars on the MAMI/DAPHNE points and 
the MAMI/TAPS results are only statistical. The systematic scale uncertainties are ~ 2% for 
LEGS, ~ 4% for MAMI/DAPHNE and ~ 6% for MAMI/TAPS. 
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FIG. 6. The energy dependence of the differential cross section at 6^+ = 85° for ^(7, TT~^)n in 
comparison with our energy dependent multipole analysis (full line) and the VPI solution SP97k 
(dotted line). For the LEGS points the statistical error and part of the systematic error are 
combined to the net uncertainty bars [32|. The error bars on the MAMI points and the Bonn 
results [53 1 are only statistical. The systematic scale uncertainties are ~ 2% for LEGS, ~ 4% for 
MAMI and ~ 6% for Bonn. 
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FIG. 7. The differential cross section da/dQ and the photon asymmetry S for p(7,p)vr'^ (left) 
and p{'y,7r~^)n (right) at = 320MeV. The MAMI data are compared to the results of LEGS 
and our energy dependent multipole analysis (full line). For the LEGS points the statistical error 



and part of the systematic error are combined to the net uncertainty bars [32|. The error bars on 



the MAMI points |26| are only statistical. The systematic scale uncertainties for da/dQ are ~ 2% 
for LEGS and ~ 4% for MAMI and ~ 2% for the photon asymmetry of MAMI. 
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FIG. 8. The difference of the hnear polarization cross sections, {da±/d^ — da\\/dCl)/{2sin^{6)) 



T,da / dQsin ^9 for p{^-,p)t^^ and p{^, 7r^)n. 
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FIG. 9. The contribution of higher partial waves {It^ > 2) to the differential cross section for 




(MeV) 



E^ (MeV) 



7P piT^ at 6*^ = 0°, 6*^ = 90° and 6^ = 180°. The solid line is the ratio (da'P /daf""^^) of the 
differential cross section for only s and p waves to the differential cross section including higher 
partial wave contributions (Born contribution for > 4). The curves are obtained from SAID |39] 
solution SM95. 
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FIG. 10. The sensitivity of differential cross sections da/d^ and photon asymmetries S to the 
Ei^ multipole and higher partial waves {l^ > 2) forp(7,p)7r° and p{j,TT+p)n at = 340MeV. The 
dashed hnes are obtained with s and p waves only. The dotted lines excludes the amplitude 
and the full lines show the results with all partial waves. 
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FIG. 11. The isospin 1 = 1/2 and / = 3/2 components of the s- and p-wave multipoles for 



the proton |^6[ in compared to the energy dependent result of the fixed-t dispersion analysis |27|. 
The plotted quantities ReA^y. /cos(5y_ = A^f^e"**^'* are equal to the absolute value | A^f± | up to 
a sign, see eq. (4). 
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FIG. 12. The real and imaginary isospin / = 3/2 components of the Mi+ and Eij^ multipoles. 
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The solid dots show the result from the energy independent fit only to the MAMI data (da/dQ 
and S) [^. The solid line and open circles show the energy dependent and energy independent 
results from the fixed-t dispersion analysis [|27[. 
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FIG. 13. The energy dependence of the ratio E-^j^ • The sohd dots (energy independent) 

and the sohd hne (energy dependent) present the result from the fixed-t dispersion multipole 



analysis ||2^. In addition the energy dependence of R = C||/(12A||) is shown as open squares. 
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